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ABSTRACT 

We study the fields of spherically symmetric thin shell sources, a.k.a. Dyson spheres, in a fully 
nonlinear covariant theory of gravity with the simplest galileon field. We integrate exactly all 
the field equations once, reducing them to first order nonlinear equations. For the simplest 
galileon, static solutions come on six distinct branches. On one, a Dyson sphere surrounds 
itself with a galileon hair, which far away looks like a hair of any Brans-Dicke field. The hair 
changes below the Vainshtein scale, where the extra galileon terms dominate the minimal 
gradients of the field. Their hair looks more like a fuzz, because the galileon terms are 
suppressed by the derivative of the volume determinant. It shuts off the 'hair bunching' over 
the 'angular' 2-sphere. Hence the fuzz remains dilute even close to the source. This is really 
why the Vainshtein's suppression of the modifications of gravity works close to the source. 
On the other five branches, the static solutions are all singular far from the source, and 
shuttered off from asymptotic infinity. One of them, however, is really the self-accelerating 
branch, and the singularity is removed by turning on time dependence. We give examples of 
regulated solutions, where the Dyson sphere explodes outward, and its self-accelerating side is 
nonsingular. These constructions may open channels for nonperturbative transitions between 
branches, which need to be addressed further to determine phenomenological viability of 
multi-branch gravities. 
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1 Introduction 

Recently there have been numerous attempts to deform General Relativity with the addition 
of extra operators which reside in the infrared (for a review, see [1]). The main purpose was 
to see if such deformations may yield novel avenues to address the cosmological constant 
problem. In the process, it was realized that such deformations generically bring in new 
degrees of freedom in the gravitational sector, by breaking residual gauge symmetries of 
General Relativity. The new theories, which are not trivial rearrangements of a scalar-tensor 
gravity in some cumbersome gauge, display a phenomenon of strong coupling, which at the 
classical level leads to the breakdown of the perturbative description of the extra modes at 
very large distances, on backgrounds with heavy sources. A prototypical example of such 
behavior was already encountered in Fierz-Pauli massive gravity [2], where the extra mode(s) 
do not decouple even when the graviton mass is taken to zero [3] , and where nonlinear effects 
become big at distances that can be much larger than the Schwarzschild radius of the source 
[4]. The latter observation, made by Arkady Vainshtein [4], prompted him to suggest that 
the nonlinear effects might conspire and tame the extra modes in the graviton spectrum at 
distances shorter than the scale where linearized perturbation theory breaks down. Similar 
behavior has been argued [5] to also arise in other models of modified gravity such as DGP 
[6], galileon theory [7], and the new model of massive gravity recently proposed by de Rham 
and Gabadadze [8] (see, for example, [9, 10, 11]). 

However a complete demonstration of how and why this effect works is lacking. There 
were attempts to integrate background equations in purported nonlinear extensions of Fierz- 
Pauli theory [12, 13], but they have conflicting outcomes. In [12], it was noted that the 
solutions which are flat at infinity do not connect onto the smooth interiors near sources, 
but become singular in between the desired asymptotic regions. More recently, in [13], nu- 
merical exploration suggested otherwise, and those authors claimed that the solutions which 
extrapolate smoothly between mass sources and asymptotic infinity do exist. Nevertheless, 
since this work was based on what is really not a complete and fully covariant theory, one 
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may worry if the additional terms may have been omitted which could change the outcome 
of the numerics. Indeed, even those authors warn of the sensitivity of the perturbative ap- 
proaches on nonlinear effects. Moreover, even if that is not so, a more explicit and direct 
understanding of this phenomenon is warranted, as much of the work on nontrivial deforma- 
tions of General Relativity at some level utilizes the so-called "Vainshtein mechanism" . One 
would therefore wish to have a foolproof demonstration of operability of the said mechanism, 
in a way that leaves little doubt as to what may be going on. 

That is the purpose of this work. Here we shall explore the simplest theory which fea- 
tures the Vainshtein mechanism, the fully nonlinear covariant cubic galileon, with conformal 
couplings to matter [14]. To simplify the analysis, we will focus on the fields of spherically 
symmetric thin shell sources, or Dyson spheres [15]. Thus the effects of matter sources on 
the fields will be modeled by the boundary conditions which the fields satisfy across matter 
shells. This simplification does not reduce the generality of our results, but it does enable us 
to integrate exactly all the field equations once in the static case. Hence the resulting field 
equations governing static solutions reduce to first order nonlinear equations. For the cubic 
galileon, static solutions come on six distinct branches. 

On one, a Dyson sphere surrounds itself with galileon hair, which far away looks like a 
hair of any Brans-Dicke field, with its usual long strands, ~ 1/r, governed by the standard 
4D Gauss law. The hair changes at the Vainshtein scale, from wherein the extra galileon 
terms dominate the minimal gradients of the field. In contrast to the long canonical hair, 
their hair looks more like a fuzz, because the galileon terms are suppressed by the derivative 
of the volume determinant. This shuts off the 'hair bunching' over the 'angular' 2-sphere at 
shorter distances. Hence the short fuzz is less dense even close to the source. We stress that 
this is different from screening in gauge theories, where the medium develops spontaneous 
polarization which shields the fields generated by sources. Instead, it's the nature of the 
hair which changes, while the Gauss law formula and the force felt by probes remain the 
same. The transition between the two hair regimes is controlled by the dimensional coupling 
constant of the cubic galileon term, which sets up the Vainshtein scale ry. When the mass 
scale of this coupling is chosen to be low, the distance ry is long, implying that the hair 
fraying occurs far from the source, where the hair is dilute and hence the extra force field 
from the scalar is weak. The diluted fuzz ensures that the field strength can not grow very 
large further inward. This is really why the Vainshtein's suppression of the modifications of 
gravity works close to the source, as we will discuss in detail. In turn, the price one pays is 
that the effective Lagrangian for the galileon is cut off at a low scale, albeit we confirm that 
in terrestrial conditions this cutoff is close to a millimeter if one picks the galileon couplings 
such that they could yield interesting cosmological effects at very large distances. 

On the other five branches, the smooth static Vainshtein's coiffure does not stick on. All 
the static solutions on these additional branches are singular far from the source, and closed 
off from asymptotic infinity. We show this qualitatively, but the proof is general, as we use 
the phase space analysis of the system, where the singularities manifest as the bifurcations of 
the phase space trajectories. Four branches of singular solutions are descendants from similar 
singular solutions in Brans-Dicke gravity [16, 17]. Yet, the fifth is really the self-accelerating 
galileon sector where the singularity is a manifestation of forcing the geometry to stay static. 
Static solutions in the self-accelerating sector are inconsistent, since the scalar condensate 
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far away drives time evolution. This is reminiscent to what happens, for example, in the case 
of singular global strings [18], which become regular when allowed to inflate in the core [19]. 
The difference is that while the scalar hair for global strings is fully supported by the string 
core, here the local sources merely trigger the scalar gradient. The faraway evolution is 
controlled by the playoff of the scalar terms, with only minor contributions from the sources. 
Although we do not analyze the Vainshtein phenomenon on the regulated time dependent 
self-accelerating branch, we expect that in the simple model which we analyze cannot really 
yield a healthy dynamics there, because the scalar field on that branch is a ghost. So while 
the short range hairs of this ghost may be suppressed, its long range behavior would still 
yield pathologies. On the other hand, if the ghost may be excised by further modifying 
the theory, and leaving the self-acceleration on smooth backgrounds, we believe that the 
Vainshtein effect may remain operable on such branches too. Yet, it remains to be seen if 
multi-branch theories without any ghosts anywhere actually exist. 

Motivated by the similarity with singular global strings, we construct explicit examples, 
in thin wall limit, of regular time-dependent solutions in our covariant cubic galileon theory. 
In such configurations, the Dyson sphere is really an expanding bubble wall, separating self- 
accelerating branch from a (section of) asymptotically flat branch, and vice versa. Such 
configurations can be viewed as bubbles of one vacuum formed inside another, and are 
channels for nonperturbative transitions between branches. Our bubbles indicate that the 
self-accelerating branch is unstable to decay to the asymptotically flat branch, with the 
transitions described by nucleation of bubbles with positive tension. The bubbles describing 
pockets of self-accelerating branch inside the Minkowski vacuum, however, require negative 
tension, and it remains to see if such bubbles can be consistently resolved within the energy 
spectrum of the theory itself [20, 21]. This is particularly interesting because such non- 
perturbative channels may occur in other modified gravities with strongly coupled modes 
and multiple classical branches of solutions [20] , particularly if they reduce to galileon models 
in the decoupling limit. Hence it is important to study this issue further, to determine 
phenomenological viability of such models. The galileon models offer a setup within which 
this important and interesting question may be addressed efficiently. 

2 Setup 

As we already stressed in the introduction, we will work with the cubic galileon theory with 
conformal couplings to matter, and restrict matter sources to reside on the surface of Dyson 
spheres. This is described by the following action [14], 

S = J d 4 x | (^R - ^(90) 2 - d(90) 2 V 2 0) - V^Anatter^, matter fields) | , 

(1) 

where g^ = e^^^ Mpi g^. We choose the overall sign of the leading order galileon term, 
— 1(<90) 2 , to ensure a stable ghost-free Lorentz-invariant vacuum. It then follows that the 
so-called self-accelerating vacuum is plagued by ghosts. In addition, we assume, without loss 
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of generality, that a > 1 , but allow f3 to take either sign. This action can now be split into 
the "bulk sector" , off the Dyson spheres, and the sector residing on the spheres, where the 
equations simply reproduce Israel conditions relating the geometries inside and outside the 
Dyson sphere. Away from them, the simplest approach to setting up the bulk problem for 
the static background is to work in the action on spherically symmetric backgrounds. When 
the background geometry is taken to be static and spherically symmetric, 

ds 2 = -e 2u dt 2 + e 2X dr 2 + e 2a d{l 2 , 
= <P(r), (2) 

we can dimensionally reduce the bulk action for the gravitational sector with the galileon on 
the Killing vectors d t , dg, d^. A straightforward calculation yields the following ID effective 
action 

S™ k = Jdr {(2<rV + a' 2 y +2 °- x + e» +x - I e "+ 2 -y 2 - ^(v* + 2<y')^e u+2 ^^ . (3) 

The primes denote r-derivatives. We have divided the whole action by M Pl , redefined the 
scalar by = MpiLp, the galileon coupling by a = aMpi, and dropped the irrelevant overall 
factor of angular area 4ir. In arriving at Eq. (3) we have used the fact that the cubic 
galileon term, y/^g(d<p) 2 V 2 (p is really only a product of first derivatives, up to an irrelevant 
boundary term, because on the background (2) it becomes ^e~ 2u ~ 4a (ip' 3 )' . 

Varying this action with respect to the metric functions v, A, a and the scalar <p we find 
the equations of motion in the bulk. Indeed, the full system is 

(e" +2ff -y)' + 2a (V + 2a')<p /2 e u+2 °- 3X )' = , 

y?' 2 + 4a(u' + 2 ( x')y? / V 2A = 2(2vW + a' 2 ) - 2e 2X ~ 2a , 

^ e u+2a-\y _ a ^>3 e u+2a-3Xy = 1 J + (J /2 )e ,+2a-A + 1 g ,+A 1 ^+2a-A,„/2 



3 vr 2 V 2 4 

"(// • 2rr') r -'V" 
3 



[(„> + o'y +2 °- x x - — y v +2CT - 3A y = {2v'o' + cj l2 y+ 2 °- x - V+ 2 -y 2 

3 2 

-— (v' + 2a')<p' 3 e» +2 °~ 3X , (4) 
3 

where the equations are results of varying (3) by ip, A, u, a, respectively. It is obvious that 
the former two equations are first order (after trivially integrating the first equation). The 
latter two look more formidable. However, subtracting twice the third equation from the 
fourth, we find 

[(*/ - a')e v+2 °- x }' = -e v+x . (5) 

Now, since we are at liberty to gauge fix the "lapse function" A at will by picking an 
appropriate coordinate transformation r — > /(r), we can choose it such that \ = —v. In 
this gauge Eq. (5) reduces to \{y' — a')e u+2a ^ x }' = —1, and so it is also readily integrable. It 



1 If a < we let (p — > ~(f> and (3 — > — j3 to recover a theory with a > 0. 
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suffices to keep only it and the first two of the system (4) because by the covariance of the 
theory, the four equations in (4) are not all independent. One of them is a linear combination 
of the (derivatives of the) others, as enforced locally by Bianchi identities and stress-energy 
conservation. 

Therefore in the A = — v gauge the full system in the bulk reduces to a first order system 
of three differential equations, 

e 2 "+ 2 V + 2a(u' + 2a') V l2 e Au+2 ° = Q , 

y?' 2 + 4a{u' + 2a'y z e 2v = 2(2v' a 1 + a' 2 ) - 2e- 2v ~ 2 ° , 

v' = a' - e- 2l/ - 2a (r - q) , (6) 

where Q and q are two integration constants, to be determined by the boundary conditions 
to follow. The fourth equation from (4) is redundant as noted, being a linear combination 
of these, and their first derivatives. 

On to the matter sources on Dyson spheres. As we noted, they just induce sources in 
the Israel junction conditions relating the bulk configurations inside and out. Thus we can 
incorporate them in the standard way. However, in the spherically symmetric static case, 
we can use the shortcut, by employing the Einstein-galileon field equations with a 5-function 
matter sources and using the Gaussian pillbox trick to relate the singular terms in the bulk 
equations and compute the jumps in field first derivatives across the shells, while demanding 
continuity of the metric and scalar fields. The reason is that in the spherically symmetric 
static systems, the metric and scalar (2) are already written in almost the Riemann normal 
form. Further, we only need to obtain three independent boundary conditions because the 
second equation in (6) was already of the first order in the bulk form (4), and so it does not 
yield nontrivial junction conditions. 

To proceed, we note that the shell stress-energy in our gauge is given by 

T M j,(matter) = diag(— p, 0,p,p)S(r — r ) , (7) 

where p and p are the energy density and the tangential pressure at the shell, and r the 
shell radius. We have normalized the energy scales in the matter sector on the shell relative 
to the local time coordinate t = e u °t, because this absorbs away the local redshift of the 
matter sector coming from the gravitational potential well induced by the shell's own mass. 
In other words, p and p are proper energy density and pressure. The radial pressure vanishes 
because, as we noted, the radial Einstein-galileon equation is first order in derivatives, as in 
ordinary General Relativity. The two nontrivial conditions, come from the "00" and, say, 
the "22" equations of the following, 

G\ = T» v ((p) + ^(matter) , (8) 

where we have normalized all dimensional scales to Mpi. They are given by, 

2a"e 2u = -pS(r - r ) + ^e 4i V 3 )' + regular terms , 

{y" + a")e 2u =p5(r- r ) + — e 4 *V 3 )' + regular terms . (9) 
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Meanwhile the scalar field equation, 



VV + « (2V M (VVW) - V 2 (Vv?) 2 ) + |T%(matter) = , 



(10) 



yields 



(e^+V)' + 2a ((!/' + 2a / )v? /2 e 4!/+2,T ) / + |e 2<T (2p - p)5(r - r ) = . 



(11) 



Integrating (9), (11) across the shell, between r ± e and taking e — > finally gives us the 
required junction conditions. With all the functions evaluated at r = r , they are 



Here the symbol AP designates the difference of the quantity P between the different sides 
of the Dyson sphere. We should stress that both the interior and exterior solutions must be 
appropriately gauge fixed to ensure that the Dyson sphere resides at the same value of the 
radial coordinate viewed from either side. 

In any case, we see that to determine the gravitational and galileon fields surrounding a 
Dyson sphere we need to find solutions of (6) satisfying the junction conditions (12) on the 
sphere. Let us turn to that now. 

3 Growing Hairs: Bulk Solutions 

Note that the solution that describes the interior of a regular static stable Dyson sphere 
should be unique, a spherical section of flat Minkowski space with a constant galileon tp. The 
reason is that the static galileon configurations which continuously deform the Minkowski 
vacuum do not violate null energy condition. If we had turned them on in the interior, they 
would have induced a curvature singularity in the center, as dictated by Hawking-Penrose 
theorems. The singularity could be avoided by the null energy violating, time-dependent 
galileon configurations, but such backgrounds would be plagued by ghosts if the Minkowski 
vacuum is ghost-free, as in the example that we focus on. Therefore, the interior solutions 
degenerate to empty space. 

This should not be too surprising, since after all one encounters precisely the same situa- 
tion in conventional General Relativity in the presence of normal matter sources. Moreover, 
such configurations are also suggested by black hole no-hair theorems, which require that 
black holes can't support primary scalar hairs, setting their charges to zero. Quite so, we 
note that the Schwarzschild solution is a solution of the cubic galileon theory if we set the 
integration constant Q to zero, and pick the trivial root of the first of Eq. (6). Given this, 
we expect that such standard black holes should form by the collapse of Dyson spheres, if 
it is controllably slow and if during it scalar hair is shed away by radiation. It is hard to 
imagine that an interior which contained a lot of scalar hair could shed it quickly enough 




Av' + Aa' = e' 2u p + ^e 2 * A(y?' 3 ) . 



(12) 
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to let a completely regular black hole form. This reasoning is already highly suggestive that 
Vainshtein configurations should exist somewhere in the theory once we pick empty interiors. 

Let us first consider the limiting case a — in (6). In this case the theory reduces to the 
simple scalar-tensor Brans-Dicke theory, parameterized by the single coupling to matter /3. 
The bulk equations (6) are then exactly integrable, as is well known [16, 17]. The solutions 
in our gauge A = —v can be written as [17] 



where £,5,^ are integration constants, the latter two subject to the condition 5 2 + 7 2 = 1. 
Note that for 7 7^ 0, the solution (13) in actual fact is a compact way of writing two separate 
families of solutions: 

• one family extrapolates between the curvature singularity at r = £ (which is really the 
"center" of the coordinate system, as seen from the second of Eq. (13), because there 
the area of the constant r spheres vanishes) , and Minkowski geometry at infinity; 

• another family extrapolates between two curvature singularities, r = and r = £. 

In this case the theory has the symmetry tp -h- —tp, reflected in the solutions as the invariance 
under simultaneous transformation 7 — > —7, the radial shift r — >■ r + £, and the sign flip 
£ — > —£. This symmetry is trivial on the first family of solutions at infinity, where the scalar 
ip vanishes. The solutions which are regular at infinity flow to the fixed point of the discrete 
symmetry ip -h- —(p, given by ip — 0. To tailor together the full geometry, we need to pick the 
family of solutions based on asymptotic properties far from the source, and use the junction 
conditions on the shell to uniquely specify the specific solutions within this family. The latter, 
doubly-singular, family of solutions can be thrown away in conventional gravity, be it General 
Relativity or Brans-Dicke, because one of the two singularities would be far from the source, 
and the solution would not flow to the vacuum outside of matter. Junction conditions, 
specified with usual matter that satisfies conventional energy conditions, never forbid us 
from doing this. Furthermore, quantum effects do not reintroduce the singular solutions, 
because the transitions from regular backgrounds are suppressed by large Euclidean actions 
demanded by the singular regions. However, in the case of modified theories, in principle the 
junction conditions may alter these conclusions, and one may have to pay special attention 
to this issue. 

Let us now return to (6) with a > 0. To classify the families of solutions which the 
system admits, we note that the discrete symmetry <p -B- — tp is explicitly broken by the 
simultaneous presence of quadratic and cubic galileon terms. As a result, the degeneracy of 
the solutions which belong to the same families of the quadratic theory will be lifted, and 
they will be rearranged in six new families of solutions, in general: 




(13) 
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• the family which extrapolated to regular Minkowski infinity will splinter into two new 
families. One will still behave as the parent family far from the source, with negligible 
cubic contributions, and tp flowing to zero. In the other family the cubic contributions 
will remain comparable to the quadratic ones, such that they continue to nearly cancel 
far away, leaving a small mismatch that corresponds to the nonzero value of Q, as seen 
in the first of Eq. (6). This family must have a larger tp than the previous one, which 
therefore cannot flow to zero far away. Thus this family will be singular far away; 

• the doubly singular family will now split four ways. The reason is that it was char- 
acterized by different asymptotic limits for tp at the singularities r = and r — £, 
and tp didn't vanish at either end. So, first the symmetry tp -H- —tp will be explicitly 
lifted, and furthermore the cubic terms on each individual branch will flow differently 
depending on the relative signs of the quadratic and cubic galileon contributions on the 
two branches. Hence, in the static limit we will see four different families of solutions, 
all singular at both ends, but interpolating between the two singularities in different 
ways. However, one of these families can be regulated far away if we allow it to become 
time-dependent, since it corresponds to the self-inflating galileon branch. 

This is precisely confirmed by investigating (6). From the first and third of Eq. (6), we 
can express v' and a' as functions of <p', 

3 1 q}+ 6a V ' 2 e^+ 2 - ' 

a' = - e - 2u - 2 °(r-a) + - ^ (14) 

3 1 q> 6a^' 2 e^ ' 1 ' 

and then substituting these two in the second of (6), after some straightforward rearrange- 
ments, obtain 

(Q - e ^+ 2 V) 2 6ay 4 e 4 " 
Q 2 + |(r - q) 2 + 2e 2 ^+ 2CT ~ 1 + 6aV 4 e 4 ^ ' ^ ' 

This is a sextic in <//, and for the general case of Q ^ it has six different roots. These 
roots define the six different branches of solutions of (6), which are generated by taking 
any individual root of (15), substituting tp' as a function of v, a and r on it into (14), 
and integrating these two equations. Note, that the expression for tp' doesn't need to be 
integrated any further, because the actual value of tp is not physically relevant, thanks to the 
pseudosymmetry of the theory (1) defined by the simultaneous shift of tp and the rescaling 
of dimensional parameters in the matter Lagrangian. 

We cannot solve (15) exactly, and write the specific roots in closed form. Fortunately, 
that is not necessary in order to understand how the solutions behave. We can determine the 
qualitative behavior of the solutions graphically, by considering the roots as the intersections 
of the two curves defined by the different sides of Eq. (15) at any fixed value of r, and 
consequently of v and a, and then vary the radial location and the metric functions to 
see how the intersections behave. So, the left hand side (LHS) of (15) for fixed r, is, a is a 
parabola in 0', shifted from the origin by Qe~ 2u ~ 2a . The RHS defines a curve which is a 
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quartic parabola near the origin, which changes to a slowly growing function saturating at 
unity as ip' — > ±00. This function looks like a smoothed square well, because very close to 
the origin the quartic parabola is very shallow. In a generic situation, with Q 7^ 0, the two 
curves are offset relative to one another, so that they touch the horizontal axis at different 
locations. Therefore they must intersect at several points. 

At radial distances close to a Dyson sphere, where the exterior metric coefficients are 
not too big - specifically, where the area of the angular sphere is not much larger than the 
area of the Dyson sphere itself - the two curves will intersect six times, as illustrated in Fig. 
(I). What happens is that e v = 0(1), because we are taking Dyson spheres to be larger 




Figure 1: Roots of Eq. (15) defining the branches of the cubic galileon theory. 



than their Schwarzschild radius, so that their fields are not very strong, and near them the 
gravitational redshift is small. Simultaneously e u = 0(r), and so the coefficient of <p' on the 
LHS is 0(e 2u+2a ) ~ 0(r 2 ), and the shift from the origin is 0(Q/r 2 ), being typically very 
small. On the other hand, while the quartic is very shallow near the origin, its coefficient 
is 0(a 2 ), and because the length scale yfa is very large, and so as <p' increases, the RHS 
changes sharply and quickly saturates at unity. 

Now we can follow the system as the radial coordinate is increased, moving away from 
the Dyson sphere. Notice, that in this case the behavior of the system is very sensitive to 
which root we select, because the variation of the functions v and a is controlled by (14). 
However: even without a detailed investigation, we can easily deduce that five of the six 
roots must be singular at some distance away from the Dyson sphere. Indeed, if we increase 
the functions v and a in (15), the potential well shape in Fig. (1) narrows down, while the 
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parabola moves towards the origin and becomes very steep very quickly, trying to climb off of 
the potential well. Therefore the pairs of roots 3&4 and 5&6 are moving closer together, until 
they touch and eventually disappear altogether. This means that they cannot exist after v 
and a increase beyond some critical finite value, defined by the coincidence of the two roots. 
So, in turn, v and a cannot grow indefinitely if we are on one of these four roots. Hence 
these four families of solutions cannot asymptote to flat Minkowski space, where a does grow 
beyond any bound. Furthermore, the moment when the roots osculate is a singularity of 
the system of differential equations, because the two pairs of families of solutions become 
degenerate. This means that the sphere of the radius where the degeneracy occurs must be a 
curvature singularity outside of the Dyson sphere, just like in the case of the doubly singular 
solutions in Brans-Dicke gravity (13). Note also, that we can estimate the critical distance 
where this occurs, by noting that the roots become degenerate when each side of the Eq. 
(15) is roughly 0{\). This yields, with e CT ~ r, to <// ~ 1/r from LHS and if' ~ l/^/a from 
RHS, and so r cr ~ y/a. 

On the other hand, roots 1 and 2 might allow v and a to grow indefinitely on them, as r 
grows larger. To see what happens we can try to venture very far from the Dyson sphere on 
these two roots, and check if the fields are self-consistently getting weaker. To start, at very 
large r we may try to approximate if' <C 1/ \fa and e a ~ r to the leading order, in order to 
check if these two roots are remaining close to the origin in Fig. (1) when r, z/, a are growing. 
Substituting this approximation into Eq. (15) yields 

(q _ e 2-+2<y) 2 = 4a 2 e 4l - + V 4 (l + 3e 2 ^) 

where we have used the fact that the denominator on the LHS of Eq. 15 goes like 2r 2 /3 + 
2e 2u+2<T ~ 2e 2<T /3 + 2e 2u+2a . Since to leading order if' ~ Qe- 2u ~ 2a , setting e = Q- e 2u+2a ip' 
and plugging it here, we find 

e 2 ~ 4a 2 Q 4 e- 4 "- 6CT (l + 3e 2 ") . (16) 
So the roots 1&2 of the equation (15) are given by 

^ = Q e - 2 - 2CT (l-2aQ e - 2 ^ 3CT (l + 3e^) 1/2 ) , <p' 2 = Qe~ 2v - 2 ° ' (l + 2aQe~ 2u - 3 ° ^(l + 3e 2 ^) 1/2 ) , 

(17) 

which follows because the root 1 is closer to the origin than the root 2. Now we can substitute 
these two roots into Eq. (14), finding 

v' = -\e- 2 ^{r-q)±^{l + ^) ll2 e- 2 ^ , 

a' = i e - 2 - 2 -(r-g)±^(l + 3e 2 ^) 1/2 e- 2 ^, (18) 

where the top sign corresponds to the root 1 solutions, and the bottom sign to the root 2 
ones, respectively. 

With the top sign, the function v has a fixed point at infinity, where v = 0, and the terms 
0{l/r) in the top equation cancel precisely, leaving u' = 0{\/r 2 ). Likewise, in this limit 
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a" — 1/r + 0(l/r 2 ), consistent with all of our approximations. The asymptotic behavior of 
the root 1 solutions at large distances is 



e 2 ^ 1 - 



(19) 



r 



like the standard Schwarzschild solution and its scalar-tensor generalizations (13). The root 
1 is certainly moving closer to the origin, and so it continues to exist. Therefore, the branch 
defined by the root 1 is indeed asymptotically flat, and has fields which become very weak 
from the source. 

With the bottom sign, however, we find that sufficiently far away v must stop growing 
and begin to decrease. The first of (18) dictates that v' must become negative. The key is 
that on this branch the scalar field gradient <// is slightly greater that Qe~ 2 "~ 2a , and this 
is independent of the approximations made in getting Eqs. (18). So, since near the Dyson 
sphere e u is positive, this implies it should vanish somewhere far. On the other hand, a may 
continue to grow, as controlled by the second of Eq. (18). Since <// = Qe~ 2v ~ 2a to the leading 
order, on this branch the value of ip' will grow, instead of decrease as on the asymptotically 
flat branch defined by root 1. So our approximation ip' <C l/y/a will break down, and we 
need to go back to the full Eq. (15), where we should now consider the limit where r, a 
and p> ! are growing, while v is decreasing, to determine what's really going on. In this limit, 
the coefficient of the <p' 2 on the LHS of (15) is e 2l/+2a /r 2 , and it is decreasing, while the 
location of the minimum, Qe~ 2u ~ 2a is increasing. So the parabola is getting shallower and 
rapidly moving away from the origin, farther to the right. Meanwhile, on the RHS, since e v 
is decreasing, the potential well is getting wider. This means, that the roots 2 and 3 of Fig. 
(1) are moving towards one another as r increases. At some distance they converge together, 
and become degenerate, which signals that they are singular there. So the solutions from 
the branch characterized by the root 2 will also have a naked singularity far from a Dyson 
sphere sourcing them, again at a distance r cr ~ y/a. 

It is instructive to consider various branches of the cubic galileon from the perspective of 
the decoupling limit, which is widely employed in galileon models, as well as other modified 
gravity setups. The decoupling limit is formally defined as the limit where gravity is de- 
coupled by taking Mpi — > oo, effectively fixing the background geometry, while keeping the 
scalar, its self-interactions, and any matter which sources the scalar. Since we have absorbed 
Mpi into ip and a, to enforce the decoupling limit in our equations we can merely rescale 
ip — > Zp, a — >■ ol/Z and Q — > ZQ, take the limit Z — > oo, fixing e 2u = 1 and e°" = r. The 
equations in (6) degenerate to a single scalar field equation, 



These two roots define the branches of solutions given by roots 1 and 6 of the general case. 
Note that the normal behavior familiar from existing decoupling limit considerations require 



r 2 ip' + Aarip' 2 = Q , 



(20) 



with only two roots, 




(21) 
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Q > 0. If Q is negative, then the description based on (21) breaks down at r* ~ (a\Q\y^, 
a scale which as we shall see shortly is just the Vainshtein scale of the source. 

The upper sign gives the root </?'—>■ % in the limit of large r. It corresponds to the branch 
defined by the root 1 of the full theory, which admits a flat space limit, and interpolates 
between the near region of the Dyson sphere and Minkowski space at infinity. We stress 
again, that for the root 1 solutions to be really able to extrapolate between infinity and 
the surface of a small Dyson sphere with weak fields, the discriminant of the root equation 
must be positive, and Q > 0. If this is not the case, the Dyson sphere cannot develop large 
Vainshtein regions, as we will discuss in more detail in what follows. 

The lower sign gives the root <p' — > — ^, a confining solution, which yields the branch of 
solutions defined by the root 6 of the full theory. We can see this from (15), where in the 
decoupling limit the RHS reduces just to the quartic, implying that the roots 4 and 5 cannot 
be observed in the decoupling limit at all. Since the second root of (21) is negative, it must 
be root 6. These solutions are really self-accelerating galileon backgrounds in disguise, as is 
clear from the decoupling limit 2 . This is rather telling regarding the nature of the singularity 
at large distances on this branch. While it is certainly present for static backgrounds, since 
they have a confining potential far away, its appearance is a consequence of imposing the 
static limit. Consistent with this, the estimated location of the singularity at r ~ ^Ja is 
precisely at cosmological scales where modifications of gravity alter the vacuum, and where 
y/a plays the role of the Hubble scale. When galileon is allowed to be time-dependent, the 
geometry is completely nonsingular far away and in the future. Thus time dependence is 
a kind of 'censor', removing the singularity on this branch. We will consider this in more 
detail in what follows. On the other hand, the asymptotically flat solutions are deformations 
of the flat vacuum induced by localized sources. Hence, it is those solutions which should 
be featuring the Vainshtein phenomenon, when we consider gravitational fields of Dyson 
spheres on that branch. We will examine it in detail in the next section. We also see that 
the decoupling limit analysis is blind to the existence of other, strongly gravitating solutions. 
This is not surprising, because we know that those solutions reduce to (13) in the limit a — > 0, 
and those cannot be deduced without gravity. 

Finally, we note that this analysis confirms our assertion that the interior of Dyson 
spheres must be an empty space, because all the solutions which have nontrivial scalar hair 
have short distance singularities. As those are unacceptable for regular matter distributions, 
the interior must be free of scalar hair, and hence flat. 

4 Vainshtein's Coiffure 

The branch defined by the root 1 of (15) yields solutions which are asymptotically flat. 
So, as noted above, this is the one candidate family of configurations which should feature 
Vainshtein mechanism. Let us examine them and find out exactly what the Vainshtein 

2 Note that if we tried to enforce the decoupling limit by 'blunt force trauma' approach in (15) we 
would have found a square of (20) - i.e. r 2 (p' ± 4ar(p' 2 — Q - admitting also roots 2 and 3. This is 
inconsistent, however, since these two branches do not exist outside of gravitational sector, as evidenced by 
their singularities. Technically, in obtaining (15) we have employed gravitational equations of motion, and 
encountered cancellations between powers of Mpi, that yield spurions in the 'decoupling limit' of (15). 



12 



mechanism is. In the literature it is often denoted by a term Vainshtein's screening. This is 
a misnomer, in our view. If we were to take the idea of 'screening' literally, it would imply 
some effect of 'renormalizing' the source charge similarly to what happens in electrolytes with 
Debye screening, or in quantum field theory. There, the external fields are compensated by 
the induced polarizations from real or virtual dipoles, arising in response to the external 
fields, and reducing the net charge of the source. Normally, the dipoles require negative 
charges, which is quite fine in gauge theory but possibly disastrous in gravity, where negative 
masses yield instabilities. On the other hand, gravity could be neutralized by 'free fall', but 
here what needs to be suppressed is the forces due to the extra mode(s) (whose charge is 
still the mass, so the quandary of negative masses persists). 

The resolution is, as we find it, that in fact Vainshtein effect is really due to the dynamical 
redefinition of the long range hairs sourced by a mass rather than 'conventional' screening. 
The net charge remains the same regardless of where it is measured, as controlled by the 
simple vacuum Gauss law. Instead, the field strength which it sources changed dynamically 
because of the irrelevant operators which become big closer in. Quantitative aspects of it, 
however, have been correctly anticipated in [5], which nevertheless assumed the decoupling 
limit. Here, we will in fact derive this in a self-consistent manner from the theory given 
by (6), on the branch defined by root 1 of (15). This leaves no doubt as to the existence 
of smooth interpolating solutions across the Vainshtein's radius. Further, our construction 
reveals the limitations on the types of matter that feature the Vainshtein effect, revealing 
that the Vainshtein radius does depend on the specific microscopic details of the matter 
source, as anticipated in the formulation of the so-called elephant problem [22, 1]. In this 
section, we will sketch the basic ideas which explain the Vainshtein effect, relegating the 
more precise, quantitative statements about matter featuring it to the next section. 

So, let us start by looking at the limit a = 0, which is the pure Brans-Dicke theory, 
in order to define the problem. The extra scalar couples to matter through the conformal 
combination g^ v = g^e^ . This is the cause of the problem. In perturbation theory, this 
gives the coupling 0(l)(3(pT^ M , in our normalizations. This means that a mass M surrounds 
itself with the scalar hair. By the standard 4D Gauss law, the scalar hair is 

= < 22 > 

because the 'charge' of the source is 0(1)/3M/Mpi. Consequently, the energy density carried 
by the hair outside of the source is p v ~ M 2 Pl y' 2 = 0(l)f3 2 J ^. In this conformal frame, 
the scalar energy density can't affect the gravitational potential of the source mass which is 
still given by the Newtonian expression Vn = — Gn t m . However, probes couple directly to 
the scalar hair, because they move on geodesic trajectories of the conformal metric g^ v , and 
not of g^y. The matter trajectories are defined by 

d oc^ i , doc doc „ / , 

and so the force from a spherical mass distribution which acts on non-relativistic matter 
probes is F nct = Tq . In terms of the original metric g^ u , using the standard conformal trans- 
formation formula with the conformal factor f2 2 = e^, this gives F net = Tq — |^p-- Then, 
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using T r 00 \d r h m d r V N = -G N M/r 2 and d r Q 2 /Q 2 ~ /V = 0{l)(3 2 M/{M P l 2 r 2 ) = 

0(l)/3 2 G N M/r 2 , we find that the net force is 

F net = -G N ^(1 + 0(1)(3 2 ). (24) 

If the effective coupling /3 2 is not smaller than about 5 x 10 -3 , the theory would run into 
problems with solar system bounds, coming from the Cassini satellite. 
Turn a back on. The Gauss law formula (22) changes to 

p' + 2a{v' + 2a VV = 0{l)^e~ 2 »- 2 ° , (25) 

m pi 

as per the first of Eq. (6). Here we have used Q = f3(p — 2p)r\j2 = 0(l)(3M jMp 3 , found 
after inserting the first of Eq. (12) into the first of Eq. (6). We take Q > 0, necessary to 
interpolate weak fields to short distances, as noted previously 4 . For Dyson spheres larger than 
their Schwarzschild radius, r > r s ~ G^M, on this branch we can consistently approximate 
e 2v ~ 1 — - and e CT ~ r. 

r 

So to the leading order we find precisely the decoupling limit equation for the scalar, 

, det q' (9 . . n M , . 

where we have rewritten 2a' = 2 = in a more illustrative way, as the logarithmic 

derivative of the metric determinant. This changes how the field lines distribute over a fixed 
solid angle at shorter distances. Indeed, if the configuration is asymptotically flat, at very 
large distances the Gauss law fixes the field line source to be ~ M/r 2 . There asymptotic 
flatness ascertains that the cubic terms are tiny and the hair (26) is approximately equal to 
(22). As the distance decreases, the hair grows thicker, the cubic terms catch up, and begin 
to compete with the standard ip' field strength. This occurs at p' ~ a^f^p' 2 , or therefore 

where p' ~ 0(l)/3 M f 2 ~ r/a. At shorter distances, the cubic contributions completely 

pi r 

take over. The transition scale is precisely 



rv* I ^V /S - MP\f 3 , (27) 



L PI 



i.e. it is the Vainshtein's scale 5 . Below ry this yields a different scaling for p'; an additional 
effect comes from det g' / detg ~ 2/r in Eq. (26), which dramatically slows down the bunch- 
ing of field lines at shorter distances. This term partially closes the 'angular' 2-sphere to the 
galileon field lines, which therefore have a larger spread over a given solid angle at shorter 



3 Hcrc, 0(1) term disguises the sensitivity of the hair to the type of matter distribution, i.e. the elephant 
problem, alluded to by the p — 2p <~ 1 — 2w prefactor, which we will discuss more in the next section. 
4 We will also discuss this in more detail in the next section. 

5 Actually, it depends on the type of matter source, being proportional to p — 2p <~ 1 — 2w. 
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distances. Taking Eq. (26) into account for Dyson spheres which are much larger than their 
Schwarzschild radius, at distances below r v , we find 

✓ - om {^-T - °w f — f - w 1 (^) 

yaMpjr / \ ar J r v V r 

The field lines are /ess dense at short distances than they would have been without the 
galileon contributions. Their distribution is depicted in Fig. (2). The galileon suppression 
by the higher dimension operator involving the derivative of the volume determinant makes 
them look more like a fuzz, spreading them out and preventing 'hair bunching' closer to the 
source. Further, because the length scale y/a is taken to be large, the transition scale for 




Figure 2: Vainshtein's coiffure. 



this 'fraying', or 'splintering' of the scalar hair ry is very large compared to the gravitational 
radius of the source. Hence where the transition occurs the field strength is very weak, and 
so it remains very weak all the way down to the source, thanks to the modified distribution 
of field lines (28). 

This suppresses the galileon corrections to forces acting on probes inside the Vainshtein 
radius. First, they do not gravitate very much. The scalar energy density is p ~ a^f-cp' 3 ~ 

3 

o:(3 2 ^f (^r) 5//2 . Because rs < r <C ry, p < /3 2 (^-) 1,/2 ^, this is much too small to directly 
alter the background geometry, as can be readily checked by comparing it to the RHS of the 
second of Eq. (6). Since the nontrivial gravitational corrections scale as ~ T^/r 3 , they are 
orders of magnitude larger than the galileon sources, shifting the metric by terms oc l/y/r, 
which dilute more slowly than 1/r, but are tiny to start with. Second, and as before, since 
matter moves in the background given by g^ v = e^g^, the effective force a particle following 
the geodesies (23) is F nct = Tq — The conformal factor derivative now is still oc cp', 
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but the scalar gradients are different. From (28) we find d r Vl 2 /Vl 2 



/V = G(l)/3 2 ^ 




Therefore the net force acting on nonrelativistic probes at distances below the Vainshtein's 
scale is 



Integrating this once, the net nonrelativistic gravitational potential a particle inside the 
Vainshtein radius is subject to is 



The correction is AV/Vn = /3 2 (^) 3 ^ 2 , and, leaving aside the galileon coupling (5 and the 
type of matter correction oc 1 — 2w, this is precisely what was found in DGP model, in 
[23]. This confirms that the cubic galileon model can be viewed as a decoupling limit of the 
simplest 5D codimension-1 setup in DGP [6], on normal branch where ghosts do not roam 
[20]. The consequences of the correction are small, as noted in [23, 24], with the dominant 
effect being the precession of bound orbits, as per the (violation of) the Bertrand's theorem, 
which implies that any deviations of central forces that differ from ~ 1/r 2 or ~ r yield orbits 
that are not closed. This has been argued to be relatively small, but only by an order of 
magnitude or so below the current experimental sensitivity of lunar laser ranging [25]. In 
fact, for the case of galileon, the coupling parameter f3 gives a better chance for observability, 
and the simple cubic galileon may already be constrained by the lunar laser ranging, as we 
shall discuss below. 

The analysis above explains the argument that the Vainshtein effect is a consequence of 
long range hair fraying rather than 'conventional' screening as in gauge theories. The total 
charge, defined by the usual AD Gauss law never changes, but the field lines which it sources 
are redistributed between the extra field strengths included in the irrelevant operators that 
take over at sub- Vainshtein distances. 

In closing our expose of the Vainshtein effect, let us note its implications for the pertur- 
bative description of the scalar fluctuations around the background. The key is to determine 
the full quantum cutoff, which as noted in [26] is the scale beyond which the gravitational ef- 
fects are probing the trans-Planckian distances. Since the scalar perturbation theory naively 
breaks down at the scale where the cubic scalar terms in (1) take over from the quadratic 
ones. This occurs at 1/a 1 / 3 [26], or equivalently at the distance scale r cuto fr — {^pict) 1 ^ 
(which is parametrically the Vainshtein radius of a Planckian black hole). If we were to take 
y/a to be Hubble's distance, y/a ~ Hq 1 ~ 10 61 £ P i, so as to ensure that it is only cosmo- 
logical solutions and little else that are affected by the modification of gravity, this would 
have yielded the r cuto fr — 10 41 £p/, or therefore r cuto fr — 10 3 km. Such a low scale is the main 
source of problems in DGP, as noted by [26]. 

However, what really breaks down at this scale is the linearized description of scalar 
perturbations around the the far-field vacuum, at super- Vainshtein scale distances. Below 
the Vainshtein scale one needs to account for the background effects of the scalar hair which 
gives larger contribution to the scalar propagator than the vacuum terms. This renormalizes 
the scalar perturbation kinetic term by a factor Z ~ atp" ~ a<p' /r, which is large for r < ry, 




(29) 




(30) 
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and pushes down the galileon self-coupling, given by the cubic scalar term in Eq. (1). After 
canonically normalizing the fields, this renormalizes the canonical coupling a — > a/Z 3 ^ 2 . 
[27, 28]. The scalar cutoff is the energy scale below which the renormalized cubic is smaller 
than the quadratic, and so it is given by A cuto fr — (Z 3 ! 2 ja) x l 3 ~ Z l l 2 {M P ij a) 1 ! 3 . So the 
short distance cutoff r cutoff ~ l/A cutoff is r cutoff ~ (i P ia) 1 / 3 /Z 1 / 2 . Using (27) and (28), the 

/ \ 3/2 

suppression factor is Z ~ ( , where r v and r are the Vainshtein radius of the source 
of the galileon hair, and r the distance where the hair strength is evaluated, respectively. 
Thus, at a distance r from the center of a source of mass M, the effective cutoff is 

,,,, / r \3/4 / a \i/i2 / r 3 \i/4 
r cu „M <PlQ )V3(_) =( w ) ( m ) . (31) 

In terrestrial conditions, with y/a ~ Hq 1 ~ 10 61 £«, this gives 

lcm 

r cutoff — 1^11/4 ' ' ' 

of order the Earth's Schwarzschild radius. This simple analysis generalizes the revised DGP 
cutoff given in [27]. While this cutoff is very low, a price to pay for having a large length scale 
y/a ~ Hq 1 , which yields large Vainshtein's radii and efficient suppression of extra forces at 
distances shorter than them, while having a chance at yielding branches where cosmological 
effects of the galileon may play a role, it can be pushed below the current experimental 
bounds, at around (9(100) microns in the gravitational sector [29], by increasing the coupling 
to matter \/3\. Indeed, it is clear from (29) that we can increase \/3\ by a couple of orders 
of magnitude without affecting the force significantly in the laboratory conditions, because 
the Vainshtein radius can be kept large and the hair near the source thin. Note that such 
increase of the extra force would not affect weak equivalence principle bounds, since it is 
independent of the type of matter probes. On the other hand, it may also happen that even 
though the galileon sector ceases to be predictive below a centimeter, whatever resolves it 
cannot produce sufficiently dramatic effects to influence any terrestrial experiments at sub- 
millimeter scale, thanks to the huge suppression generated by the large Vainshtein radius of 
the Earth, (cm/r^®) 3 / 2 ~ 10~ 28 , as argued by [27]. Thus one could take a view that while 
galileon models can't make any observable predictions for tabletop searches for submillimeter 
deviations from Newton's gravity, they may at least be safe from any bounds that follow 
from tabletop experiments. 

So, ignoring theoretical prejudices, this shows that in galileon models with more tunable 
parameters controlling the interactions may at least give a chance for searches of their sig- 
natures by correlating cosmological and astrophysical searches. As we noted above, lunar 
laser ranging is already constraining cubic galileon theory, because the precession of lunar 
perigee here may be faster than in DGP. The precession rate is controlled by the relative 
correction to the Newtonian potential [25], and in our case the advance per orbit is 

H = 3 i^ = m 3/2 H DG P, (33) 

4 V N 

which means that increasing the matter coupling (3 by only a factor of a few, while keeping 
other parameters unchanged, brings the cubic galileon effects to observable levels. This 
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yields a constraint f3 < O(10) if ^fa ~ H 1 . It would be interesting to see a more thorough 
investigation of possible galileon signatures. 



5 Combing the Coiffure: Boundary Conditions 

So far we have been discussing the solutions "in the bulk", branch by branch, in order to 
gain insight in the qualitative behavior of geometries corrected by galileon hairs. Here we 
complete the discussion by showing how to construct the realistic fields of Dyson spheres by 
using the junction conditions (12) to connect the interior and the exterior into the complete 
configuration. After we take the interior of the Dyson sphere to be a ball of Minkowski space, 
we must pick the same coordinate gauge as the exterior metric, and further ensure that in 
this gauge the Dyson sphere sits at the same coordinate location r as measured from the 
outside. 

So, starting with ds 2 _ = —di 2 + df 2 + f 2 dil 2 , and recalling that in our gauge the exte- 
rior metric is ds 2 + = —e 2u dt 2 + e~ 2u dr 2 + e 2cr dQ2, we see that the appropriate coordinate 
transformation in the interior is t = e u °t, f = e~ u °[r + (e u ° — l)r ], where r = e ao , and the 
subscript zero implies that the functions bearing it are evaluated at r = r . Physically, this 
coordinate transformation simply means that we have to synchronize the clocks just inside 
and outside of the Dyson sphere, accounting for the gravitational redshift at the sphere rel- 
ative to an observer at infinity. Following this, however, we must shift the radial coordinate 
to compensate for the Lorentz contraction of the radial coordinate in our gauge, induced by 
the redshift. This gives 

ds 2 _ = -e 2u °dt 2 + e- 2uo dr 2 + e~ 2v ° (r + {e v ° - l)r ) 2 dtt 2 (34) 

for the interior metric. Therefore, on the shell we have 



o U 



e 



e = e- u , e a = r , v'_ = (p'_ = , a'_ = . (35) 

Substituting this in the junction conditions (12), and subtracting the second from the third 
of (12), yields the derivatives just outside, 

V?' + + 2ae 2 »« (v> + + 2<) ^ = \e~ 2v \p - 2p) , 

+ r 2 + 3 V+1 

< = e- 2 -^ + fe 2 V+ 3 . (36) 

Combining the first of these equations with the first of (6), and the difference of the third 
and second with the third of (6) yields, respectively, 

Q=\{p- IvYl , Q=(P + PYl + (1 - e<> . (37) 

The first of these two relations has been used in the previous section already. 
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The second of Eq. (36) is the condition which fixes the radius of the Dyson sphere, once 
its energy density and pressure, and the branch of the exterior solutions are given. This 
equation may have several roots, or none, depending on the branch and p and p. The second 
equation of (6), when evaluated at the Dyson sphere, 

tpl + Aa{v' + + 2a' + )<p , y° = 2(2u' + a' + + < 2 ) - . (38) 

does not yield further significant constraints on the parameters of the solution. When ip' = 0, 
its RHS has a nontrivial eigensolution e 2v = 1 — £/r, e a = r, on which it vanishes. This is just 
the Schwarzschild solution, and we could always cut it at an arbitrary radius after specifying 
the right energy density and pressure at the shell. These many degenerate solutions mean 
that this equation merely controls the deviation of v' away from the Schwarzschild form as 
governed by the scalar sources, rather than fixing the radius of the Dyson sphere. In the 
weak field limit, we can therefore ignore it as redundant. 

So now in principle one can work out the exterior derivatives for any of the bulk branches 
and find the full exterior geometry. Without galileon modifications, we normally keep the 
solutions which asymptote to flat space as the physically relevant ones, and throw away the 
ones singular far away as having unacceptable boundary behavior. Hence the most pressing 
check is to see whether the galileon terms alter this in some way. 

On the one hand, as we have already indicated previously, when a > 0, the sources with 
P(p — 2p) < prevent weak fields from setting up inside the Vainshtein sphere. To see it, it 
suffices to consider the decoupling limit, where the root 1 branch is determined by the upper 
sign solution of Eq. (21). Using the first of (37) yields 




When a/3(p — 2p) < the discriminant vanishes precisely at r ~ (a|/3(p — 2p)|rg) 1//3 ~ r v . 
If we define the shell equation of state, p = wp, this gives ry = (a\(3(l — 2w)\pr^) 1 ^ . The 
weak field description below the Vainshtein radius fails, and if the discriminant vanished 
for any finite value of r, by continuity we would expect a branch merging where uniqueness 
of solutions breaks down, which signals a singularity. So a Dyson sphere would have to be 
placed farther out from any such surface, when the stress energy on it satisfies the condition 
a/3(l - 2w) < 0. 

On the other hand, when the sources have stress energy satisfying p < —2p, the shells 
cannot be static. To see this, we can treat the system (36) perturbatively, ignoring the 
scalar field at first. In that case, the regular exterior solution reduces to Schwarzschild, 
e 2v = 1 — ijr and e a = r. The boundary conditions (36) then give 

Note, that the first formula reproduces the correct expression for the active gravitational 
mass of a source [30]. Clearly, if we take p + 2p < 0, £ must be negative by this equation. 
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But then the second equation can only have a negative solution for energy density. For 
sources obeying null energy condition, p + p > 0, this requires p > — p = \p\, yielding 
p + 2p > 0, which contradicts the starting assumption. Thus, barring direct null energy 
violations in the matter sector, as we claimed we must have p > — 2p, yielding £ > 0, for a 
static shell to have any energy on it at all. In fact, it is easy to understand this intuitively. 
The point is, that if we were able to make a static Dyson sphere with sources satisfying 
p + 2p < 0, we would have obtained an exterior solution which would look like Schwarzschild 
with negative mass, and without singularities in the core. That would have yielded many 
problems in the theory, starting with instabilities and so on [31]. Hence, such solutions 
should not be possible, and indeed the junction conditions exclude them. 

When we turn on the scalar, the conclusions do not change for asymptotically flat solu- 
tions, r > £, if we ignore the galileon terms, because the junction conditions for the metric 
remain qualitatively the same after we substitute (13) into (36). Once we turn on the galileon 
corrections, the extra terms on the RHS of the second of Eq. (36) could, for <p' > 0, trans- 
late into additional positive contributions on the RHS of the second of Eq. (40). One may 
worry that such corrections could yield static Dyson spheres with positive energy density, 
but obeying p + 2p < 0. However, this is very difficult and technically requires huge Dyson 
spheres, that may be very unstable and subject to large IR corrections. We can see this as 
follows. Suppose that a is very large, so that Vainshtein effect may be operable. It makes 
the extra galileon corrections tiny, ip' < 1/ y/a. So, any solution with p > and p + 2p < 
must have at most a tiny energy density p ~ 1/ ^fa. When the matter does not violate the 
null energy condition, the inequalities p + 2p < and p + p > further limit the pressure, 
p ~ 0{p). So this gives £ ~ 0(rl/ y/a), and that, from the second of Eq. (40) with the 
galileon corrections included, would yield tq > 0(y/a), making the Dyson sphere larger 
than the scale where gravity is modified! In other words, such spheres would have to be 
cosmologically large. So in a cosmological setting, the cosmological effects, and causality in 
the very least would obstruct their existence 6 . 

As a bottomline, we see that the 'conventional' Vainshtein effect may be operable for 
sources on the Dyson sphere which satisfy the following energy conditions 

af3(p-2p)>0, p + 2p>0, p>0. (41) 

For these sources, static Dyson spheres may be built, and their external fields are to leading 
order conforming to general relativistic expressions at distances below the Vainshtein radius. 
At larger distances, the fields reduce to the linearized solutions with scalar fields (13), behav- 
ing as in a scalar-tensor theory. So, for sufficiently large a, giving a large Vainshtein radius, 
the theory may be in accord with observational bounds, and the scalar forces may indeed 
be negligible, if we take solutions on the branch of root 1. In order to recover the results of 
conventional General Relativity, specifically for sources which behave as nearly pressureless 
dust, from this equation we see that a (3 should be positive. This means, that the sources 

Nonetheless, it is curious that at least in theory one might imagine constructing such huge negative mass 
static sources, even on a branch of the theory where there is no ghost. In fact, one may use the existence 
of such solutions to place a bound on y/a. If one wishes to leave no chance for any late formation of such 
large spheres with repulsive Newtonian potentials, and prevent any chance of them wreaking havoc in a late 
universe, one concludes that yfa > H^ 1 . 
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with more negative pressure cannot benefit from the protection of the Vainshtein effect, and 
the fields of extra modes will not stay small near such sources. Hence arguments that the 
cosmological volumes of matter must automatically be subject to the description appropriate 
to the 'interior' of the Vainstein sphere [5] may be incompatible with the operability of the 
Vainshtein effect at shorter scales if the net effective pressure inside the cosmological region 
is sufficiently negative to reproduce cosmic acceleration. 

Note, that the Vainshtein radius depends not only on the total mass of the source, but also 
on the type of matter which is sourcing the fields. For matter satisfying (41) the Vainshtein 
radius is given more precisely by 



which shows that the suppression of extra forces only works when pressures are not too 
large. As the pressure increases, the Vainshtein radius shrinks, while simultaneously the 
scalar hair, and the force it generates, became smaller. Once the scalar gradients become 
negative, however, the suppression of the hair is impossible for asymptotically flat solutions, 
since the two terms in the galileon conservation law get locked together, and so the hair does 
not splinter enough. As we noted, this supports the idea that the Vainshtein effect must 
be sensitive on the details of the source which generates the fields rather than just its total 
mass, as suggested by the elephant problem. 



The discussion of the previous section shows that Dyson spheres which carry stress energy 
with < p < —2p won't remain static. That is as it should be, since generically such 
shells are spherical domain walls separating different vacua of the matter sector, and will 
eventually expand out, generically becoming null at late times, with p — > —p. In conventional 
General Relativity, these domain walls separate bubbles of lower vacuum energy from the 
parent vacuum with higher energy density, and are widely used for cosmological applications 
(see, e.g. [33]). In modified gravity models with multiple branches, a new possibility arises, 
where bubbles may catalyze transitions between branches. Here we shall outline this for the 
cubic galileon model. Nevertheless we think this possibility is fairly generic, and needs to be 
understood in more detail. 

In what follows, we will make several simplifying assumptions, since we only aim to 
demonstrate that branch transitions are possible in principle. For this reason, we will work 
with the assumption that the vacuum energy on each side of the bubble is the same, and 
that the energy difference between the sides is provided by the galileon gradients. Further, 
we will work perturbatively around the decoupling limit, which simplifies the calculation, 
and at least allows us to consistently describe the asymptotically flat and self-accelerating 
branches. The latter assumption in itself is not too restrictive for large bubbles, but the 
former one is more special, and in general one would like to extend the analysis by adding 
generic nonvanishing cosmological terms off the shell. We hope to return to this elsewhere. 

Bubble mediated branch changes have been considered in the cubic galileon in the decou- 
pling limit previously, in [21]. These authors found an instanton describing the nucleation 




(42) 



6 Hirsute Bubbles: Branch Changelings 
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of a tensionless bubble separating two branches, but then they argued that such solutions 
were merely a quirk of the thin wall limit, and that tunnelling between branches was in fact 
prohibited. Their point was that the vanishing effective tension was a sign of exact cancel- 
lation of positive and negative energy densities, where the latter had to be put in by hand 
to mediate the transition, and was otherwise extraneous to the theory. We remain dubious 
about this interpretation. The self-accelerating branch contains a perturbative ghost [20], 
and generically one expects such branches to be always present in multi-branch models [32]. 
So such ghosts may provide the requisite negative energy in the core of the wall, with a 
perfectly physical matter sector contributing positive energy. Obviously, this implies that 
the self-accelerating side of the bubble must be modified relative to the geometry used in 
[21], but that may be necessary anyway once the ghost is accounted for. Motivated by this 
reasoning, we will reconsider the issue of branch changing bubbles here, and, in the leading 
order gravity approximation, provide solutions that describe nucleation of hirsute bubbles 
with large and positive tension, which circumvent the no-go claims of [21]. 

Let us now show how to construct such branch changelings in the theory (1). We take 
the parent state to be a vacuum, either on the asymptotically flat branch or on the self- 
accelerating branch, and the descendant to be the state with a bubble of the other vacuum 
surrounded by the parent vacuum. For simplicity, we will set vacuum energy contributions 
from the matter sector to be zero, and work only with the self-accelerating contributions 
from the galileon. This suffices for our purposes here. To describe a bubble, we take an 
SO (3, 1) symmetric configuration 

ds 2 = g^dxPdx" = e 2n{0 (dC? + CHjd^dx^ , <p = (p(() , (43) 

where r ) i jdx % dx^ is the metric on three-dimensional de Sitter space with unit curvature. We 
take the bubble wall to lie at a fixed "radius", ( = ( w > 0, such that the field profiles are 
given by 

lU K«*(C) C>C» ' ^ g lv^(C) C>C« ' 1 ] 

The interior and exterior are taken to lie on opposite branches. This ansatz is just a deforma- 
tion of the standard VIS geometry used to describe tensional domain walls in flat space [34], 
which was utilized in the construction of inflating branes in braneworld models [35]. The 
'appearance' of staticity in (43) is therefore just a mirage, coming from the codimension- 
1 structure of the bubble wall. In the reference frame of an inertial observer at infinity, 
however, one can easily ascertain that the bubble is expanding. 

To find the field profiles in (43), we can again work in the action as before. It is then 
convenient to change the gauge (i.e. 'un-fix' the gauge) in (43) by defining new 'radial' 
coordinate, by dz = e~ 9 d(/(, where 6 is an arbitrary function, and redefine the field 7r to 
7r = 7i + In (. Transforming the ansatz (43) to this gauge, and substituting it into the the 
action (1), after factoring out the common powers of Mpi, we get 



_ (Te 3(if+j8 V >/2)| 



*(C»)}> (45) 
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where we have performed some integration by parts, taking care of any relevant boundary 
terms. The second line of (45) is the bubble wall tension term at z — z(( w ), with physical 
tension a = M Pl a, given our normalizations. To obtain the field equations governing the 
configuration (43), we now vary this action with respect to <p, 6 and it, and after obtaining 
the field equations, gauge-fixing to 9 = 0, tt = n + In (, eliminating z for (, and integrating 
the galileon ip field equation once, we obtain, respectively, 



0, 



e 2 V + 6V 2 Q + tt'J = p , 

6C 3 e 2 V + C 4 3e 2 V 2 - ^e 2 > /2 - Qa V ' 3 ^ + vr'^j 

^(^^^(CV^ + ^^e.' 2 -^ 2 ), (46) 

off the domain wall, where, as before, Q is an integration constant which follows from shift 
symmetry <p —> (p + constant. Primes designate ^-derivatives. On the wall, the the junction 
conditions relating the sections of (43) inside and out are given by the continuity conditions 
An = Aip = 0, and the analogue of the Israel equations 
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where as before A represents the difference of an appropriate quantity between the outside 
and the inside of the wall. These are easily obtained from (46) with the wall contributions 
included as 5-function terms following from varying (45), and using the Gaussian pillbox 
integration, as before. Further, as we discussed in previous sections, on the wall we can 
always absorb the value of f into the definition of the matter coupling , and so we can 
always set ip — on the wall. 

Exact solutions of this system of equations are not known to date. However we can 
employ an expansion in 1/Mpi, and substitute into (46), (47) the perturbative ansatz 

and recall that a = Mpia. To leading order we recover the decoupling limit in which 
Mpi — > oo with a and a/Mpi held fixed. The higher order gravity corrections are then 
obtained iteratively, order-by-order in Mp]. After the calculation is done, we can absorb 
the one extra power of Mp} to render the leading order solutions dimensionless, as per our 
conventions. The key is to find the solutions which are smooth on either side of the bubble. 
The reason is that such solutions have smooth Euclidean continuation, and so one would 
expect to have an instanton describing their nucleation which has a finite action. This is 
made easier when one turns on matter vacuum energies inside and outside the bubble, and 
allows them to dominate asymptotically far from the bubble. In that case, one will find 
instantons which are merely "squashed" standard tunneling solutions in General Relativity, 
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looking like creased spheres. In our case, what suffices is to demonstrate that even without 
matter vacuum energies nonsingular configurations exist. We will not go after determining 
their action here, although that task warrants further investigation. 

The configuration describing a bubble of self-accelerating branch in the vacuum of the 
asymptotically flat branch in the interior has the fields 

«' U <-0 = ^(d-?) + O(M- F ?), d -f = -l + 0(M- p }), (49) 
while outside they are 

^\0 = O{M P ?l ^ = -Y2^ak 3+0iMp ' ) - (50) 

In deriving these solutions, we have imposed regularity at £ = 0, and the asymptotic bound- 
ary condition 7r out — >■ 0. This suffices since away from the wall we can always shift tp by a 
constant, as we explained above. The galileon charges are given by 

Q m = 0(M P ?), Q out = -^^ 2 +0{M- p f), (51) 



whereas the brane tension is 



'•3 



°=-w^ + ° {M " ] - (52) 

For any finite size bubble, the tension is negative. This suggests that to nucleate the bubbles 
of the self-accelerating branch in the asymptotically flat vacuum of the normal branch may 
be difficult, suggesting the need for ghost matter fields to account for the negative tension. 
Specifically, if we were to try to associate the tension only to the wall-localized matter, we'd 
need to put negative energy modes there. On the other hand, since there is a ghost already 
in the theory, on the self-accelerating branch, one might imagine that in some resolved, thick 
wall description of the transition, the kink of such a ghost induced by the bubble wall of 
normal matter, might provide the requisite negative energy. Indeed, the scaling of the field 
gradients in the solution with the size of the bubble are curiously in accord with such a 
possibility. However, without a more detailed demonstration that such kinks may actually 
arise, which is beyond the scope of the present work, it is impossible to pass the immediate 
verdict on the transitions in this direction. Furthermore, when the matter vacuum energies 
are turned on, this problem will get even more complicated thanks to the difference of the 
vacuum energy on different sides of the wall to the wall's tension. Those contributions 
may even make the overall tension positive, and allow for transitions to the self-accelerating 
branch even without negative energies on the wall, or near it. This issue deserves further 
scrutiny, and we hope to return to it in the future. 

Describing the reverse transitions is simpler. The configuration describing bubbles of the 
flat vacuum expanding in the self-accelerating background are given by the field solutions in 
the interior, 

tt^C) = 0{M~ P f) , -|- = 0{M P f) , (53) 
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which are surrounded by a self-accelerating exterior 7 



- - (4) + oiM " >• ^ = ~i- + ° (M " ' ' (54) 

Again, we have imposed regularity at £ = 0, and picked 7r(0) = by scaling £. The galileon 
charges are now given by 

Q in = 0{M~ P f), Q out = + 0{M~ P f) , (55) 

and the brane tension is 

This time the brane tension is positive for ( w > 0. Therefore transitions from the self- 
accelerating branch to the normal branch do not require any exotic types of matter to 
catalyze them. They are mediated by the coupling between the unstable gravitational sector 
and the ordinary matter fields whose excitations ultimately populate the bubble wall. We 
expect that this is a non-perturbative manifestation of the perturbative ghost instability 
on the self-accelerating branch. Note, of course, that at zeroth order in our expansion, we 
recover the result of [21] in the decoupling limit, where the tension is found to be zero. Here 
we see that the leading order gravity correction to that is one of positive tension. 

To sum up this section, we see that there exist channels for nonperturbative branch 
changing. A branch change can be accomplished locally and in a causal manner, by a 
nucleation of a nonsingular "hairy" bubble, whose wall kicks the galileon field and while 
sourcing its hair, it pushes it do a different branch from the parent configuration. Bubble 
walls which demarcate self-accelerating vacua in the interior from the external asymptotically 
flat vacuum require a negative tension, and further analysis is needed to determine if they 
can naturally arise in the theory. The bubbles which describe sections of an asymptotically 
flat vacuum in the interior from the external self-accelerating vacuum however have walls 
with positive tension, and so may be constructed using the standard matter sources. It 
also remains to see what happens when we allow for nonzero matter contributions to the 
vacuum energy inside and out of the bubble, and whether these branch changing bubbles 
really provide the means for exploring the full structure of the landscape of modified gravity, 
a view argued for in [36]. 



7 Summary 

In this work, we have been able to give for the first time the complete proof of the Vainshtein 
mechanism, which is widely assumed to operate in modified gravity models with strongly 

7 To see explicitly why this is the self-accelerating solution, note that the physical metric, as seen by 
matter, is g^ v — e l3<p g tll , which to leading order in our expansion goes like g^ ~ ^1 — r]^ v in the 

exterior region, where ( 2 — rj^x^x" . For a(3 > 0, this is just de Sitter space written as a local perturbation 
about Minkowski, along the lines discussed in section 1.1 of [7]. 
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coupled additional modes in the gravitational spectrum. Working with a simple cubic galileon 
theory and thin shell spherically symmetric static sources, a.k.a. Dyson spheres, we find the 
first integrals to all the field equations, reducing them to a more tractable first order system. 
This allows one to classify all the solutions, and specifically show that for a range of matter 
sources on the Dyson spheres, there exist asymptotically flat solutions with scalar hair. The 
scalar hair scaling with distance changes from the standard AD form, ~ 1/r dictated by 
the usual 4D Gauss law, to another scaling, controlled by the higher-dimension operator in 
the galileon sector, at a very large distance, the Vainshtein radius. This is the key to how 
the Vainshtein mechanism manages to hide the extra forces induced by the additional light 
modes in the theory from probe particles moving inside the Vainshtein radius. The hair 
changes by getting fuzzy below the Vainshtein radius, and so its field lines do not bunch up 
as densely as in the usual 4D case, remaining more dilute. Because this transition occurs 
far, where the fields are weak, they remain weak all the way very close to the source. Since 
the force acting on a probe is still oc <fi', it is therefore also very weak. 

For a(3 > 0, this Vainshtein coiffure, of fuzzy hair below the Vainshtein radius, sets up 
for matter sources on the shell whose energy density and pressure obey p > 2p > —p and 
p > 0. This includes much of the usual matter, such as radiation and non-relativistic mass 
distributions. However, when pressure is large, p > p/2, as would be the case if the matter is 
field gradients confined to the shell, the hair would not splinter efficiently, and the Vainshtein 
effect would not be operable. Similarly, for a(3 < 0, the Vainshtein coiffure is only styled for 
matter sources on the shell satisfying < p < 2p, and p + 2p < 0. It is immediately evident 
that this rules out shells composed of ordinary pressureless matter, or even a combination 
of matter and radiation whose effective equation of state lies somewhere in between w = 
and w — 1/2. All of this points to the general issue of sensitivity of the Vainshtein effect 
on the details of the source distribution, and we indeed find signs of that. Specifically, we 
show that the Vainshtein radius depends not only on the total mass of the source, but also 
its properties, as encoded by the pressure. 

We expect that the above considerations can be generalized beyond our Dyson spheres to 
more complicated matter configurations. The condition p + 2p > is just the strong energy 
condition, while a(3(p — 2p) > is equivalent to requiring the trace of energy momentum 
to take a particular sign (depending on the sign of a/3), a/3T < 0. The latter condition is 
particularly important as it controls whether or not the Vainshtein mechanism is realized 
(a(3T < 0) or if the solution is expected to run into a singularity (a(3T > 0). Note that in 
the marginal case of T = the galileon is not sourced at any scale, and it is subtle to decide 
if the Vainshtein effect is present at all. On the one hand, the absence of any galileon effects 
may suggest that for this case the Vainshtein mechanism is extremely efficient. On the other 
hand, as we see from Eq. (42) the Vainshtein radius vanishes for such sources, implying that 
the Vainshtein effect is maximally inefficient. The point is, that for such sources there is 
nothing to suppress in the first place, since they are bald to start with - the absence of any 
hair makes their fields automatically the same as in conventional General Relativity. Thus, 
in reality, for this special limit, the Vainshtein effect is irrelevant. 

In addition, we find that generically the theory admits more solutions, which in the static 
limit are all singular far away from the Dyson sphere. In many cases, one expects that those 
solutions cannot be regulated, and so may have to be thrown out of the theory, in a way 
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similar to what happens to their relatives in conventional General Relativity 8 . However, one 
of the singular solutions can be readily regulated by turning on time dependence, since it 
is really the self-accelerated, ghost infested, branch familiar in galileon models. We explore 
what happens if we try to combine such time dependent solutions with shells laden with 
sources whose pressure and energy density obey p + 2p < 0, which corresponds to the 
scalings of shell-bound 'dark energies'. In this case, we find regular solutions, where the 
shells expand separating regions of different branch vacua. These configurations point that 
the self-accelerating branch may decay to the asymptotically flat branch. The case where 
the bubbles describe pockets of self-accelerating branch inside the asymptotically Minkowski 
vacuum are more complicated because they require negative tension, and further work is 
needed to determine what really transpires. Generally, our constructions point to channels 
for nonperturbative transitions between branches, which need to be addressed further to 
determine phenomenological viability of multi-branch gravities. 
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